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The higher Sylvester waves are discussed. Techniques used involve 
finite difference operators. For example, using Herschel's theorem, el- 
egant expressions for Euler's rational functions and the Todd operator 
are found. Derivative expansions are also rapidly treated by the same 
method. 

A general form for the wave is obtained using multiplicative series, 
and comments made on its further reduction. As is known, Dedekind 
sums arise in the case of coprime components and it is pointed out 
that Brioschi had this result, but not the terminology, very early on. 
His proof is repeated. 

Adding a set of ones to the components of the denumerant corresponds 
to a succession of (discrete) smoothings and is a Cesaro sum. Using 
a spectral vocabulary, I take the opportunity to exhibit the finite dif- 
ference counterparts of some continuous, distributional properties of 
Riesz typical means. 
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1. Introduction 



This paper is a continuation of an earlier one, [1], concerned with Sylvester 
waves, Ehrhart polynomials and degeneracies in spectral problems.^ In that work I 
gave explicit formulae for the first two waves, resurrecting Sylvester's expressions. 
Here I consider the higher waves but not so far as to produce similarly clear-cut 
results. However, during the analysis, several amusing pieces of information and 
technique arose and will be exposed here. 

A classic number problem is that of restricted partitions. Given a set of non- 
negative integers d^, (i = 1, 2, . . . , d), in how many ways can a given non-negative 
integer, say /, be expressed as a (non-negative) integer linear combination of the d^l 
I call the dj, the 'components' and write this combination 

I = m^d^ = m.d . (1) 

The number of ways is referred to as the denumerant of /, following Sylvester, who 
was one of the first to study this problem in any detail. Euler gives the number in 
terms of a generating function 

/=0 ' i=l 

where the coefficient of the is the denumerant. (I use the notation of [2] to avoid 
typographical confusion.) 

It is frequently possible, in discussions involving the combination (1), to work 
throughout, more or less, with this generating function. For example, in spectral 
problems, where the denumerant might come up as a degeneracy, one could set a = 
e~* and interpret the left-hand side as a cylinder heat-kernel associated, perhaps, 
with the square root of a Laplace-like operator, after some slight adjustments. A 
Mellin transform would then yield the corresponding (^-function which, in this case, 
would be a Barnes ^-function, 



c(s,a|d) = y ^^-^ = -^ r dte-^ 



j dt{-ty-^ 



-at 



(3) 



Jl ' ' aa-e-^^*) 



^ The appearance of the present work has been delayed by factors beyond my control. 
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(On this definition, the are not restricted to be positive integers.) 

I will not pursue this global quantity further now because there is interest in 
studying the denumerant on its own, as attested to by the long history. 

Sylvester proves the theorem that the denumerant can be expressed as a finite 
sum of quantities periodic in /, called waves, each of which is associated with a root 
of unity. That associated with the root 1 is not periodic (has period infinity) and 
is a normal polynomial, W^, in I. Hence one can set. 



where U is the properly periodic part and is a sum of waves, generally written 



The whole denumerant is a quasipolynomial in I i.e. a polynomial with periodic 
coefficients, e.g. [3]. A simple proof of this was given by Wright, [4]. 

As a function of the augmented variable, I = l + ^i/'^^ ^^e denumerant I; /d; 
satisfies a parity reciprocity under I — > —I, [5] . The choice between I and I comes 
up later. 

The second wave is that associated with the root —1, and has also been 
separated in (4) because Sylvester has given an explicit form for Wi and a more 
or less explicit form for A simpler expression for W2, like that for W^, was 
presented in our earlier work, [1], and I wish here to treat the other waves by a 
similar technique. 

The evaluation of the denumerant is a linear Diophantine question related to 
the Frobenius problem and to the Ehrhart polynomial and has thus been the subject 
of extensive analysis. I refer to two books, [6], [7], for some history and motivation. 
However, Sylvester's specific formula is not referred to very often. The papers by 
Rubinstein and Fel, [8], and Rubinstein, [9], describe his method, but the detailed 
techniques are different. A discussion of Sylvester's basic method can be found in 
the book by Netto, [10] published first in 1902, although he does not give Sylvester's 
final form for W^, which seems to be largely ignored, apart from the 1909 paper of 
Glaisher, [11], who extends the working to the other waves. The present paper could 
be looked upon as an independent, partial recasting of Glaisher's computations in 
as small a compass as I could manage, plus some comments about the literature.^ 

There is some very interesting recent work on computing these partition num- 



d; 




(4) 



q>2 



bers [12], [13]. 



Unfortunately, I have not been able to see some relevant early Italian work by Trudi. 
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2. Sylvester's waves 



Sylvester's theorem leads to a prescription for the waves which is the fol- 
lowing. Write out all the factors of the components, d^. Let there be fj, such and 
call a typical one, q. Separate the components into two groups - those divisible by 
q, call them ctj, z = 1, 2, . . . o; (a is the 'frequency') and those that are not, say Pj, 
j = 1,2, . . ./3 with a + /3 = 11. Sylvester, [14], then says that. 



which stands for the coefficient of 1/t in the indicated generating function. The 
sum is over all prime qth roots of unity, p, i.e. p'^ = 1, p = e^'^^^/'^, (p, q) = 1. The 
variable I is the augmented one, I = / + ^ Q;j/2 + XI = ^ + S 

The necessary inverse power (s) of t are provided only by the group. There- 
fore, expanding the exponential, the coefficient of I /n! is 



p't 

CO 



T«,- ^n.(l-p-'^^)''''""""'^.^sinha,t/2 



(6) 



n«i V n,(l - P"''^ """" V sinha,V2 Y sin/3^.(7rp/g - it/2) 



So far as the exponent of p is concerned, it is best to retain the integer /, as Netto 
[10] suggests, also + p'?"^ + ... + 1 = 0. 

As in my earlier work, I invoke some of Euler's old products. Specifically, after 
taking logs 



smh2; ^ V n^TT^J tt^ 277^ Stt^ ^ ^ 

n=l ^ ' 
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and 



log 



sin 77(2; + a) 



sm wa 



-f^log (1 - ^) + - f^log fl + ^) 



- lim (C(l, a) - C(l, 1 - a))^ + - (C(2, a) + C(2, 1 - a))z^ 

\ (C(3, a) - C(3, 1 - a))^^ + ^ (C(4, a) + C(4, 1 - a))z^- 
1 o d II .(f 




(8) 



which is really only a check because it follows directly by school calculus. 

Reversing the argument, the step from the second to the third line, which is 
the standard reflection formulae of the Hurwitz (^-function, (or, equivalently of the 
polygamma function). 



can be derived. For this nexus of notions consult Hoffman, [15,16], who has a 
different organisation, and a little history. 



I also note that one could expand in terms of the Euler rational functions, e.g. 
Carlitz [17], [8], defined by. 



As a practical arithmetical means of calculation in any specific numerical case when 
the number of components and the frequencies are not too large, this is probably 
as good a route as any. It was used by Cayley in his treatment of partitions, [18]. 
He gives two (related) relevant expansions, based on Herschel's theorem, e.g. [19], 
[20], [21]. Firstly, for small denumerants. 



C(n,a) + (-l)\(n,l-a) 



(-1)"-^ c^^-i 
(n- 1)! rfa"-i 



TT cot 7ra , 




(9) 



n 




(10) 



1 



l-c(l + A) 
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which gives a direct expression for the Euler functions (9) in finite difference terms. 
The quantity acting on 0-^ is a finite polynomial in the difference operator, A, 
because A"' 0"^ = if n > m and, generally, 

A" 0- = (-1)- Ei-l^k-Q =^'1^1' (11) 
fe=o ^ 

in terms of Stirling numbers of the second kind, a very old result. Thus I find the 
elegant expression, ^ 

H^{X) = -(l-j^A^ 0", (12) 

which I do not take any further here but note that it is equivalent to a result of 
Vandiver, [22]. See later. 

Moreover, a form for the Todd operator, Todd(c, df^), of Brion and Vergne, [23], 

is 

a _ log(l + A) ^_„,^^ ^^3^ 



l-ce-^h c(l + A)-l 
so that 

Todd(c, a,) m = _\ fih - 0) . (14) 

For complicated denumerants, the logarithm is best expanded using (10) and 

io«T^ = /**(i-T3^) 



l-c(l + A) 



These particular expressions lead to the denumerant expressed in terms of the 
unaugmented variable, I, and this is what Cayley produces after some calculation. 
Interesting work on relieving some of this effort is given by Munagi, [13], and Sills 
and Zeilberger, [12]. 

To encounter the (preferred) augmented variable, I, a more symmetrical ex- 
pansion is required, corresponding to the cotangent in (8). The relevant expansion 
follows as 



ce-* c(l + A) 



e-0•^ (15) 



1 - ce-* 1 - c(l + A) 
Set, now, c = so that the left-hand side becomes 

1 - e-«-* 
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and the derivatives with respect to t, at t = 0, which are exphcit in (15), are just the 
derivatives of e~*/(l — e~'*) = (coths/2 — l)/2 with respect to s. This, incidentally, 
provides an elegant difference proof of these higher derivatives. 



1 dC^ f \ dC^ 1 1 + A 

coth5/2-l) = ^ i_^ = (_l)n^L±ii^O" 

2 ds" V / ds-^e^ -I ^ ^ - (1 + A) 



= -(-ir log (i-;i%i)o"^' 

n+l^ 1 



= (-1)" y - ^ -t: A'^ 0"+^ 



(16) 

using (11) to give an explicit expression. The passage from the first to the second 
line is eff'ected by the identity due to Herschel (Boole, [24], Ch.II, §13), 

(^(A) 0"+^ = (1 + A)f (A)O'^, 

whose essential content is the recurrence relation for the Stirling numbers. 

The expansion, (16) and its equivalents, surface from time to time. Agoh and 
Dilcher, [25], prove it by induction and use it, in several papers, to derive various 
Bernoulli number identities and Sterling convolutions. 

It is equivalent to a formula in Adamchik, [26] . See also Cvijovic, [27] , Knopf, 
[28], Hoffman, [15], Boyadzhiev, [29]. The technical relation used in these works is. 



d\^ „, , „ , , f ?i 1 I. d^ 

k=l 



which Knopf, attributes, in essence, to Scherk in 1824 who, so it seems to me, 
was addressing the same question solved by Herschel in 1816, arriving at the same 
result ^. Equation (17) is contained in equation (4) in Herschel, [19], and explicitly 
exhibited as Exercise 6, in Boole, [20], p. 26, viz. 



dn 1 d'^ 

dx^ k\ dx^ ^ ' 

k=l 



^ Gould gives further interesting detail and earlier history in [30] . 
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For amusement, I interject a proof of this. From Herschel's theorem one has, 



/(e*+«)=/(e^(l + A))e°-*, (19) 
and note the oft used device. 



the left-hand side of which can be picked out of (19) as the coefficient of f^/nl. The 
expansion of f{x{l + A)) in powers of A, 

/(a.(l + A)) =5;-A'=0"x^ £_/(,), 

fe=0 

after retaining powers of A no greater than n yields (18). (The first term also goes 
away because 0"' = 0.) 

These considerations also show the equivalence of (12) with a formula in Van- 
diver, [22], who also employs the operator d/dx. 

After this digression, 1 return to the computation of the denumerant, i.e. (6), 
where one has to deal with the products. The first product has been encountered 
in [1] and so I turn attention to the second one, to which (8) applies after the 
identifications, a — = /3jp/q and z = Zj = —iPjt/2. I remark that Tra^ is half the 
argument of p^i , p^i = ex.p{2ni a^) 

Denoting by Q{z,a) the argument of the logarithm in (8), the corresponding 
multiplicative sequence follows on first constructing the product, 

ft It^ It^ 1 1 1 \ 

n,n,...=exp(^--E, + -^S,--^-S3+-^3fS4 + ...j (20) 

where the S„ are defined by (8). 

S„ = z*" J cot 7ra 

= 2y(-l)";5" 0"-S n>l /o1^ 

' l-p-^.(l + A) 

l + p^^ 



after employing (16) or its equivalent. 
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In terms of Stirling numbers, cf (16), (21) reads, 

3„=2(-2.,-5:,^^|:.!{'-l}(^)', „>1. ,22) 

Things can be expressed equivalently in terms of known Euler functions. Van- 
diver has shown that, generaUy, 

dx-^\l-x) 1-x ' 
in terms of the functions defined by (9) a few examples being Hq — 1 and, 

H,(.) = -±-^, H.(^) = j^, //3W = 1J-^. (23) 

Setting X — one has that, 

(^)"(i„ot.a-l) = ^if„(/) 

Then, 

^ 2(27r)-^^^; i?,_i(p^O , n>l, (24) 

where -R^(x) are Euler polynomials, the numerators of (23). This is just (21). A 
small simplification occurs if the inversion property. 



is used. Then, 

= 2(27r)-^^/3- _ \., i?n-i(p"^0 , n>l, (25) 

which can be obtained directly. 

Any systematic analysis of Sylvester's theorem is bound to lead to similar 
quantities and Glaisher, in his extensive treatment, encounters similar polynomials. 
See [11], especially §§79,80,93-100. He also employs the difference operator A. 
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Equation (25) yields an explicit, but unsimplified, formula for the product (20), 
occurring in (6), which now has to be combined with the product over the coming 
via (7) as 

-exp(-^.^ + ^.^-^.^ + ... 

where is the sum of the nth powers of the set. ^ 
Combining the two series gives, 

f 11 
Q,iQ,2 ■ ■ ■ Q1Q2 • • • = exp I — K^t + 2^2^"^ ~ Z^^^^ ~'~ 

where ^ 

_ 1 „ 

'^2 — 1^^2~ '^2 



l_ ^„ (26) 



23 2! "3 
_ 1 „ 

^^4 — '^4 + ''■4 

etc. 

The final step expands the exponential as a power series 

9.^9.2 . . . Q1Q2 ... = I- <di{n-^)t + 02(1^1, ^2)^^ - 03(«i> «2> K3) + • • • 

where the 0^([av]j.) are all the homogeneous products of the quantities of which the 
would be sums of powers and are classic functions of the k^, (see [5], [1]). Some 
examples are 

1 1 
0o = l> 0i = «i> ©2 = 2('^i +'^2) > ©3 = -(k? + + 2^3) . (27) 

Leaving the expressions as they are, returning to (6) for the polynomial coeffi- 
cient, the polynomial for the wave takes the finite form 



^ The relation between the constants, r, here and those, r, in [1] and [5], is = 2t„. This 
reflects the fact that the series now contains all powers of t. 
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A trivial check is provided by W2, for which p = — 1 and so, since the are 
all odd, the preliminary factor in (28) is just 1/(2^ f| a). All the S„ are zero and 
the result, using (26), reduces to that in [1,11]. 

In this case the roots of unity dependence is trivial, which is not true for the 
general wave and this is the remaining computational issue. A strategic decision 
has to be taken concerning the 'final' form for the periodic polynomial. Sylvester 
in [2] writes it in terms of elementary denumerants, (I ± r); /q;, which first entails 
a reduction into a polynomial in the prime roots, p. This reduction also occurs in 
the calculations of Cayley, [31] p. 50, who expresses the final answer, equivalently, in 
terms of prime circulators, as does Glaisher, whose calculation of W^, [11] §§81-87, 
which he takes to the third term in (28) involving the square of Sj^, brings out the 
attendant complications. Andrews, [32], rewrites things in terms of the greatest 
integer function, which makes the integrality more obvious, it being obscured in 
the other formulations. In, e.g. Beck and Robins, [6], the roots of unity expressions 
are not taken further but are mostly left, and analysed, as (generalised) Dedekind 
sums.^ 



3. The simplest case of coprime components. Dedekind sums 

There is no need for complicated bookkeeping when the components, dj, are 
mutually prime. There is a subset of waves belonging to each factor on the denom- 
inator, separately. 

The expression for was given by Sylvester, [5], and was rederived in [11,1] 
so I consider the higher waves, VF^, where q divides the typical component d^, 
i = l,...,(i, also treated in [11]. Rather than write out this case specially, it is 
easier to refer to the general form (5), set a = 1, = (ij. There is only one term 
on the denominator that is proportional to t so that t can be set equal to in the 
remaining terms. This easily gives, 

_ 1 1 ^ p-' 

by setting p — )■ p~^. The are the remaining components and the p are the non- 
trivial qf-prime roots of unity, of which there are q — 1, ii q happens to be prime 

^ Beck and Robins note that Dedekind sums occur, implicitly, in the work of Sylvester and 
explicitly in Israilov, [33] in 1979. 
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when = e^'^'^/^ and also, in this case, 




Pm \ ^ Pm 




(30) 



These VF^ are just the Fourier-Dedekind sums studied in [6], and earlier refer- 
ences therein. The actual definition there used is. 




so that it is only notation to write the wave, (30), as. 



The special case expression, (29), occurs in the early short note by Brioschi, 
[34], equn.(8), and so one should add his name to those who had already encountered 
Dedekind sums. He also gives a contour proof of the general Sylvester theorem and 
I repeat it in the Appendix for the coprime case. 

Leaving the waves as Dedekind sums, (29) cannot be considered as a complete 
determination and some energy has to be expended in computing the sums. This 
is the purpose of the manipulations of Cayley, Glaisher and Sylvester, amongst 
others. Brioschi, [34], works out the particular example d = (2,3,5). Sylvester, 
[2], does (1, 2, 3), (1, 4, 7) and (1, 3, 5), for those who like variety. Remarks on the 
comput ability of the Dedekind sums are also made in [6] and cf Exercise 8.2. 

I consider the factor 11^(1 ~ P~^^)~^ where, to remind, p is a primitive qih 
root of unity and the 13 ^ are the components of the denumerant not divisible by 5, 
^ in number. It can be advantageous, as suggested by Cayley, [31], and done by 
Glaisher, [11], to reduce the mod q. Then 



q-l 



Y[{^-p-^^)-'=Y[i^-p-n 



'm 



j m=l 



where the non-negative integers depend on the set of the 13 ^ and, in general, 
have no pattern. 
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4. Roots of unity and prime circulators 



No point would be served by continuing the calculation along Glaisher's lines, 
which rapidly becomes unwieldy. Something more universal and automatic is re- 
quired. Unfortunately I cannot provide it here. 

For example, for W^, Glaisher reduces, mod q, the powers, 13 j of p before 
performing any combinations. In a general method this would be premature and 
involve unnecessary extra labour. 

Looking at the general structure (28), with (27), (26) and (21) or (25), it is 
sufficient to write each term in the summation over j in (25) as a polynomial in p 
of order q — 1 (using p'^ — 1) and then reduce any products of these (coming from 
the powers in (26)) to similar polynomials. Ultimately, the are then also such 
polynomials and one then converts Ok/YljO- ~ P~^^) a like polynomial, which 
can then be converted into Cayley's prime circulators, if desired. 

The basic algebraic problem is to reduce the ratio of polynomials of a primitive 
qth. root of unity, p^, to a similar polynomial. 



cf Battaglini (1857). Unfortunately I have not been able to see this work. It is 
mentioned with a few details in Dickson, [35] p. 121. 

At this rather unsatisfactory point I leave this detailed algebraic aspect. 

5. Denumerants and Cesaro sums 

For a given set of components, the numerical computation of a denumerant 
can proceed in several ways. The classic special case when di — 1, d2 — 2 etc. , is 
treated at length in Gupta [36], where some history is also given. The expressions 
derived by Cayley (which are essentially the same as Sylvester's as enlarged by 
Glaisher, [11]) are considered. Large I and largish d are discussed and a typical 
example is detailed. Cancellations and factorisations occur. 

Euler computed many values using recursion, but for smallish I and d it is, 
perhaps, just as easy to employ a convolution-iteration technique which yields an 
expression in terms of simple denumerants only.^ 

^ The convolution corresponds to the product form of the generating function, (2), and the 
computation to the time-honoured one of expanding each factor and collecting terms. 




.q-1 



q-1 
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The last component can be separated using the convolution 

^ ~ ^0; ^^-^^ 

di, ■ ■ ■ , d^; ^^^^ d^; d^, d2, ■ ■ ■ , d^_i; 
which is a special case of the more general, 



• • • i^d'i dj,...,d^; di,d2-, ■ ■ ■ ■,dj_i, 

Equation (31) can be iterated to the intermediate form, 



di,d2,...,dj^] ^ (^di ^ ^0 '^d-l' / =0 

1 

(33) 

and, completely, down to 

d,,d2,...,ds (^^ ds ^^^^ •••^^^ d^; cZi;- 

The simple denumerant, /g;, is 1 if g divides I and zero otherwise. I refer to 
this as Herschel's function as it is just the average of the Ith powers of all of the 
gth roots of unity, which Herschel used when introducing his circulating functions, 
[37]. In terms of the fractional part and the Kronecker delta, 

- - ^{l/Q},0 ' 

y ) 

and (34) is very easily programmed but, being recursive, is not very efficient. How- 
ever it does present the denumerant as an obvious integer.^ 

From the spectral aspect, I; /g;, Z = 0, 1, ... is the Laplacian (Neumann) degen- 
eracy on the 2g-divided circle (or, equivalently, on a 7r/g interval) (eigenfunctions, 
cos{lq9)). This can be envisaged as a one-dimensional lune. 

The convolution, (31), can be thought of as the addition of another component. 
If the component 1 is added to the simple denumerant, I; /g;, one obtains I; /l,g; 
which is the degeneracy on the two-dimensional lune, [l/q\ + 1. 



^ By this ancient brute force method, I computed 100; /I, 2, 3, 4, 5; as 46262 in 2 minutes using 
DERIVE and an Athlon 11x4 610e processor. 
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This can be repeated to give the degeneracy on the d-lune and I pursue this 
particular process, but in the context of the intermediate form (33) which yields, 

2 — r-i — : = E Z -E , . . (36) 

which I think of as a succession of smoothings of the summand or as a nested series 
of accumulated degeneracies. 

All summations can be performed except the last, and one finds 

I; (36) 

(^1, • • • , Ij-j+i; ;/=o ^ I — I' J di, . . . , dj_i; 
which is, to check, the same as (36), i.e. , 

di, ■ ■ ■ , dj-i, di, d2, ■ ■ ■ , dj_i; 

in view of the classic value of the unit denumerant. 



(37) 



A - fd-l + l' 
as follows, for example, from the expansion of the generating function, 

y — <y' = ^-,- (38) 



This denumerant gives the degeneracy on the d-hemisphere. 

The generating function version of (37) is the rather trivial splitting, 

1 1 1 



(1 - (7)^-^+1(1 - (T^l) . . . (1 - (7^^-l) (1 - (7)»-J + l (1 - (T'^l) ... (1 - a'^i-^) ' 

(39) 

It will be recognised that (36) is a Cesaro sum originally introduced to deal 
with divergent series. The standard situation (see e.g. Hobson, [38], Bromwich, [39], 
Knopp, [40]) is that for an infinite sequence, ^j,, one defines the finite sums 



4"=E(^+:-'').. 



(40) 
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and investigates the rth 'arithmetic mean' 

as n — > oo. In this expression, the index, r, although it originates as the number of 
smoothing summations, can assume any value, real or complex. Negative integers, 
however, are usually excluded (but see later). 

In terms of generating functions, (40) translates into 

oo ^ oo 

of which (39) is an example. 

The Cesaro mean is a discrete analogue of, and a motivation for, the more 
powerful 'typical mean' of Riesz, which introduces a handy continuous variable into 
the analysis, e.g. Hardy and Riesz, [41], [38]. I give only the briefest details in a 
more general setting. 

Typically, in an eigenproblem, u would be an eigenlevel label, and (7^ its de- 
generacy. A function A^(A) of the continuous variable. A, is defined as a counting 
function, encoding the spectrum A^^, 

NiX)= E 9.. 

In the special situation above, u would be /, (7^ the denumerant, /d; and the 
eigenvalue, A^, a function of /, typically a + /, where a is a constant. 
Defining the (first) accumulated degeneracy by, 

V'=Q 

or, 

one has the connection, 

N{\) = G,_^, A,_i < A < A, 

= + -gj, = Gj,_^_i — -Qj, , A = Aj^ , 
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and, of course, is the first Cesaro sum (40). 

Unfike the accumulated degeneracies, A'"(A) depends on the actual form of the 
eigenvalues, \. 

The Cesaro generating function, (41) in terms of G, is 

oo ^ oc 

^^^'^^'^ = 7r^5:^^^^- (42) 



For the Riesz mean, this iterated summation of is replaced by an (r + l)-fold 
iterated integration of N{X) which equals the Cauchy convolution, (see Knopp, [40]), 



employed specifically as a smoothing of N by Fedosov, [42], see Baltes and Hilf, 
[43], Bahan and Bloch [44]. 

Introducing the distribution (see Gel'fand and Shilov, [45] §5.5), 



the convolution is neatly written. 



N = AT * $ , 



at the same time extending r + 1 to an arbitrary variable, a, by continuation. 

Comparing this convolution with the Cesaro discrete version, (40) or (36), leads 
to the analogy 

+ + (/). (43) 

V / ; T{l + l)T{a) ^ ' 

For typographical reasons, the left-hand side is often denoted by A^~^. 

It is interesting to enlarge upon this analogy and to exhibit some continuous 
relations together with their finite difference counterparts. 

The distribution ^> has a number of basic properties, [45]. 

(i) The concentration of ^ai^) on the positive x axis. 

(ii) The convolution, 

(iii) The Laplace transform, 

dx^Jx)e-*'' = — , t>0. 



/ 

^0 



t 
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(iv) The singularities, 

Property (i) corresponds to the statement that = if Z is a negative 

integer (because of the poles in the T{1 + l)-function) unless a is a positive integer, 
when it vanishes at a finite number of negative integers, a a negative integer comes 
under property (iv). 

The convolution property (ii) is a simple consequence of the definition and 
analytic continuation. The corresponding discrete equation viz, the well known 
binomial identity, 

I 

follows likewise, most easily from factorisation of the generating function. This is 
just Cauchy's product. 

Property (iii) translates into the generating function definition, 

oo ^ 

Y.Ar^a'= a = e-^^>0. (44) 

/=o \ ^) 

Property (iv) is the most interesting one as it concerns the case when a is a 

negative integer which is usually excluded {e.g. Chapman, [46]). Looking at (43), 

the quantity of interest is, 

V{l-a) _ I r(fc + l) _ 

This also follows quickly from the binomial expansion of (44) and one could use 
generating functions systematically, e.g. Jordan, [47], to re-express the following 
remarks. ^ 

I recall now the expression for the /cth difference, e.g. Boole, [20], applied to 
the Kronecker delta, 5/ , considered as a function of I {cf Tauber and Dean, [50], 
Traub, [49] ^^). The expressions are, 

k 



m=0 ^ ^ 



I'-l-k 



s'v' = E(-ir(:)^4., = (-iy-'t!. 

m=0 ^ 



k 

I' -I 

(46) 



^ The isomorphic algebraic scheme of generahsed sequences, e.g. [48], [49], could also be employed. 
This forms an operator calculus in the field of finite differences. 
There appears to be an overall sign error in equations (3.6) and (3.7) of [50]. 
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where A is the forward difference and A the backwards one. Comparing with (45) 
the conclusion is that, 

a—^ — k 

which is the finite difference counterpart of property (iv), as I wished to show. 

Unlike the continuous version, the higher difference A 5;*^ is not concentrated 
at the origin, / = 0, but extends to the right for a 'distance' /c, emphasising the 
non-locality, or fuzziness, in its construction. As a simple feature, I remark that, 
Y^i^Q A ^i'^ — which just states that the strength of the (oscillating) 'curve' of 
a higher derivative is zero. The corresponding forwards quantity, A 5^°, extends a 
distance k to the left of the origin. 

As usual, employing forwards or backwards differences involves a manifest loss 
of symmetry which can be restored by shifting the origin of the graphs to their mid 
points. Analytically this is accomplished by the translation operator, conventionally 
written E='=-'^/^, and corresponds to using the central difference, 5, so that S'^Sf is 
the closer analogue of the Dirac derivative 5^''^{x).^^ 

Pictures of some discontinuous approximations to 5-function derivatives ob- 
tained by central differences of the (piecewise) continuous step function can be seen 
in van der Pol and Bremmer, [51], p. 83.^^ 

6. Discussion 

I have given several disparate pieces of analysis based on the explicit calcula- 
tion of a denumerant, or restricted partition. On the way a neat finite difference 
expression, (12), for Euler's functions was found and also one for the Todd operator, 
(14). An elegant derivation of a derivative expansion is likewise given. 

An expression for the general wave is given, (28), but this is not necessarily 
in its final form. 

The expansion of the denumerant into waves can be sustituted into the Barnes 
(■-function, (3), although would not aid its specific computation. 

Brioschi's early 1857 calculation of the simplest case of prime components is 
resurrected and leads to Dedekind sums, which is known more recently. 

There is, of course, an intrinsic non— locality in the continuous higher derivative. 
If the blocks there depicted are squashed to their midpoints, and e = 1/2, then these graphs 
yield precisely the symmetrical constructions, S'^S?, here. 
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Appendix A. Proof of Brioschi's formula 



For completeness, I give a pedagogic derivation of equation (29) for the wave 
in the simplest case when the components are all prime. There is nothing new in 
this, cf Brioschi in 1857, [34]. 

The starting point is always Euler's generating function for the denumerant, 

- (47) 



^=0-' aa-^'^o (i-z«)a(i-^^o 

where a is one component selected from the and the rest. Notationally I set 
a = q E Z. Extracting the power z'' using residue calculus gives 

where C circles the origin. Now blow up C to wrap round the other poles, which 
all lie on the unit circle and are given by the vanishing places of the denominator 
in (47). In this sum, I consider only those arising from from the vanishing of the 
separated typical factor (1 — 2;^) at all the non-trivial qth roots of unity, pi, . . . , Pq_i- 
This particular sum equals, by definition, the wave, W^. The contour then continues 
on to infinity, which contributes nothing. 

The residue of 1/(1 — z'^) at z = pj is easily found to be —Pj/q so that the 
residue of the integrand is pj'' /qf{Pj), and I thence obtain (29). 

The multiple pole at z = 1, coming from all factors in the denominator would 
give the first (non-periodic) wave, W^, a polynomial, essentially a generalised 
Bernoulli polynomial, e.g. [9]. There are no other multiple poles, which accounts 
for the simplicity of this evaluation. ^"^ The total denumerant is then the sum 

^.=W, + Y.^,^- (49) 

If the components are only coprime, the roots p^, for each component, separate 
into the prime roots for every factor of that component so giving a subset of waves, 
one wave, of the form (29), for each divisor. The total denumerant is still given by 
(49) with the subsum. 



Because of the separation of the roots, one might say that there is no interference between the 
waves. 
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The general case, when some components possess common factors, involves 
coincident poles on the unit circle originating from different factors on the denom- 
inator. These, like the poles at 2; = 1, give rise to polynomials. Analysing and 
organising this situation leads to Sylvester's theorem, [34]. 

Ehrhart, [52], Theorem 9.2, also gives the coprime expression for the denumer- 
ant ('compteur') which he seems to derive independently. He uses it numerically to 
compute the denumerant by calculating the expression to an adequate approxima- 
tion. The bulk of the value comes from the first wave, polynomial part. 
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